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We study in this work the ground state entanglement properties of finite XX spin-1/2 chains with 
random couplings, using Jordan- Wigner transformation. We divide the system into two parts and 
study reduced density matrices (RDMs) of its subsystems. Due to the free-fermion nature of the 
problem, the RDMs take the form of that of a free fermion thermal ensemble. Finding spectrum of 
the corresponding entanglement Hamiltonian and corresponding eigenvectors, and comparing them 
with real space renormalization group (RSRG) treatment, we establish the validity of the RSRG 
approach for entanglement in the limit of strong disorder, but also find its limitations when disorder 
is weak. In the latter case our work provides a way to visualize the "effective spins" that form long 
distance singlet pairs. 



I. INTRODUCTION 

Presence of entanglement is perhaps the most funda- 
mental difference between quantum and classical physics. 
Recently entanglement has been widely used to charac- 
terize phases and phase transitions in condensed matter 
physics. Among various ways to quantify entanglement, 
the most frequently used is the entanglement entropy be- 
tween two subsystems, obtained in the following manner. 
For a quantum system in a pure state \'ip) , the density ma- 
trix p = \il)) {'ip\ contains all information of the system in 
that state. Dividing the system into two subsystems A 
and B (most often in real space, which is the focus of this 
work), we can obtain the reduced density matrix (RDM) 
of each subsystem by tracing over degrees of freedom of 
the other subsystem: p^^^ = ti'B/A{p)- The entangle- 
ment entropy is S^ — —tr{p \ogp ) (and we can obtain 
it for subsystem B too, which for a pure state S^ ~ S§). 
While this number indeed describes how much subsystem 
A and subsystem B are entangled, it was pointed out by 
Li and Haldane^' that much more information about the 
system, especially the nature of its phase, is in the spec- 
trum of the reduced density matrices p '^ . 

It seems natural to ask if even more useful informa- 
tion is available in the eigenstates of p^/^ . The answer 
is obviously yes, as they are the basis of numerical im- 
plementation of density matrix renormalization group. 
Here we ask this question from a different perspective, 
namely we would like to ask if one can extract some phys- 
ical properties of the system from these eigenstates. We 
will show, using a very simple example, that the answer 
is still yes. The model is an antiferromagnetic spin-1/2 
chain with random nearest neighbor XX couplings, which 
can be mapped onto a free fermion model with random 
hopping. Due to the Slater-determinant nature of the 
ground state. 



P^ = 



.-H^ 



(1) 



is characterized by a free-fermion entanglement Hamil- 
tonian H^. It is easy to obtain (at least numerically) 



the single particle eigenmodes (which make up, but are 
different from the multi-particle eigenstates of p ), and 
their counterparts in part B; we call such pairs entangled 
(single-particle) modes. For the present model these en- 
tangled modes represent singlet pairs formed by effective 
spins, often over long distance (see below). 

This model is a special case of random XXZ model, 
which could also be studied using the real space renor- 
malization group (RSRGj also referred to as strong dis- 
order RG in literature).^ Within this approximation, it 
was found that the long-distance, low-energy behavior of 
this entire class of models is dominated by the random 
singlet fixed point, in which singlet bonds form between 
spins in opposite sub-lattices over arbitrarily large dis- 
tances. Refael and Moore^l calculated the entanglement 
entropy of the random singlet phase using the RSRG 
method. Within this approximation scheme, entangle- 
ment between A and B comes from singlet pairs formed 
by a spin in A and another spin in i3. If it were exactly 
correct, then in the fermion language above there should 
be a zero mode associated with each of such pair, and 
the mode wave function is (ultra)localized at the sites of 
these entangled spins. In our exact diagonalization study 
of the random XX model we find, of course, this is not 
the case unless the randomness is very strong (which is 
the regime RSRG works well); in general the (entangle- 
ment) "energies" of these modes are not zero, and the 
mode wave functions do have finite extent. We do find 
that for modest disorder strength, the entanglement "en- 
ergies" of the low-lying entangled modes (corresponding 
to long bonds that cross the subsystem boundary) in- 
deed approaches zero with increasing system size and 
bond length, supporting the asymptotically free nature 
of RSRG (including in entanglement calculation). On 
the other hand the spatial extent of these modes do not 
decrease with such increase in bond length. These mode 
wave functions are understood as real space images of the 
effective spins that form singlets over long distance; the 
profiles of such effective spins cannot be obtained within 
RSRG. Thus entanglement provides a new way to probe 
such effective spins. 



In a recent work, Fagotti, Calabrese and MooreP also 
used the free fermion mapping to calculate disorder- 
averaged moments of entanglement spectrum in this 
model. Our work is complementary to Ref. T as our 
emphasis is on entangled modes of specific realizations 
of certain random distributions of couplings, not just 
disorder-averaged quantities. 

The remainder of the paper is organized as the follow- 
ing. In Sec. II we introduce our model and notions like 
entanglement energy and entangled modes. Sec. Ill is 
a brief review of RSRG from the viewpoint of entangle- 
ment. Our numerical results are presented and analyzed 
in Sec. IV. Sec. V offers a summary and some concluding 
remarks. 



Now, we divide the system into two subsystems A and 
B, often (but not necessarily) with equal number of sites. 
We know that we can write the RDM p^ (since it is a 
positive definite operator) as the exponential of a Hermi- 
tian operator, H"^, as in Eq. (nj). For the special case of 
free fermion system here, the entanglement Hamiltonian 
H is also a free fermion HamiltonianP 






(8) 



To determine eigenmodes and eigenvalues of h we fol- 
low Ref. ii by defining correlation function 



II. MODEL AND ENTANGLED MODES 

The model we work with is a one-dimensional (ID) 
spin- 1/2 antiferromagnetic XX model with N sites and 
with random nearest neighbour couplings. The Hamilto- 
nian of the system is 



AT-l 



-" ~ 7 , -JnySriSn+i -^- S^S^^^^). 



(2) 



n=l 



Note that we do not use periodic boundary condition. By 
mapping spin operators to fermion operators via Jordan- 
Wigner transformation: 

c„ = e('^^.<"4^7)s-, (3) 

where c is a fermionic operator, the Hamiltonian becomes 



Af-l 



H = ^Yl Jn{cicn+l + cl+lC„). 



(4) 



The {N — 1) J's in Eq. Q will be generated based on 
random distribution functions to be specified later. It is a 
special case of the most general free Fermi Hamiltonian: 



N N 



H-EEh 



■ ^1 C i 



(5) 



i=i j=i 



with 



a- 



J. 



>, 



(9) 



which is an Na x Na symmetric matrix. Its Na eigen- 
modes are the same as those of h , while the correspond- 
ing eigenvalues are related to each other: 



1 



l + e"" 



(10) 



where e^ is an eigenvalue of h referred to as entan- 
glement energy (not to be confused with eigenvalues of 
the original Hamiltonian hi), and Uk is the probability of 
finding a fermion in the corresponding mode. 

We may have three different types of e. First, it can 
be a very big negative number, which means the prob- 
ability of finding a fermion at mode k is almost 1; thus 
this fermion is almost exclusively localized in A. Second, 
it can be a very big positive number, which means the 
probability of finding a fermion at mode k is almost zero; 
this means that there is a fermion that is almost localized 
completely in the complementary mode (to be discussed 
later) localized in subsystem B. These modes contribute 
very little to entanglement. Third, it can be neither of 
the above (and possibly close to zero), which means there 
is a substantial probability n less than 1 that this parti- 
cle is in subsystem A. This necessarily implies that there 
is a corresponding mode in B where the same particle re- 
sides with probability n^ — 1—n. These pairs of modes, 
referred to as entangled modes, dominate entanglement. 
To find the relation between eigenvalues of these modes, 
we use 



_r(i/2)j, ifj = * + i 

'''J 1 (l/2)J,_i ifj=*-l 



and hij = otherwise. 



(6) 



1 ^ 



(7) 



n=l 



for the ground state we expect Np — N/2 (we always 
work with even N). 



l = n' 



n^ = 



I + 6" 



l + e" 



(11) 



The number of fermions is related to magnetization: gg we have e^ 



-e^ . In the following we will use these 



relations to identify such pairs of entangled modes. 

In fact the entangled modes can be glued together to 
form a set of Np orthonormal modes occupied by the 
fermions in the ground state, and the weight of these 
modes in subsystems A and B correspond to n and n^ 
respectively. One way to obtain these modes directly 



(without using the RDMs) is using Khch's methocr^ as 
outhned below. 

We divide the Hilbert space of one particle states into 
two parts corresponding to our subsystem A and B, and 
introduce corresponding projection operators Pa and Pb 
as: 



(12) 



(13) 

where |fc)'s are the occupied eigenvectors of matrix h in 
Eq. ([5]). Then diagonalize this Np x Np matrix by a 
unitary matrix C/*^ as M ~ C/^^rf*^C/^^ , we obtain these 
new modes: 



Pa 


->;N)(^I; 


Pb = 


N 

- i. 

i=NA + l 


|^)(^ 


e a 


Hermitian matrix 








Mfefc' = 


-- {PAk'l PAk) , 





10^ = -^Y.uMPA\k} 



k 

Np Na 



k^^<"- 



kW, 



fe=l 4=1 



(14) 



III. REAL SPACE RENORMALIZATION 
GROUP METHOD 



We can obtain an approximate ground state of the 
Hamiltonian of Eq. (pi) (with J's distributed with 
a distribution function) by real space renormalization 
group (RSRG) method. As described by Dasgupta^ and 
Fisher,^ for the Hamiltonian of Eq. ([2|, we first pick 
up the maximum J — Jmax- Two spins that are cou- 
pled by this Jmax will form a singlet in the ground state 
of the Hamiltonian if we ignore their couplings to other 
spins. Then we remove the two spins that are coupled by 
this coupling constant. The two other couplings Jl,Jr 
that couple this singlet pair with the rest of chain are 
treated perturbatively, resulting in a new coupling con- 
stant J « JlJr between the further neighbor spins (see 
FIG IT]). We note that this new coupling is smaller than 
Jl and J/j, resulting in a reduction of overall energy 
scale. As we repeat this procedure many times we have 
singlet bonds that form between spins that are very far 
apart (thus forming a random singlet state) . In the mean- 
time the distribution of J's broadens (with probability of 
finding smaller effective bonds increases) , and approaches 
a power-law distribution function: 



10^ = 



VT^i^ 



WkPslk) 



Nf n 



VT^di"^ ^ 



(:^/fcC/.fe|z), 



fe=li=AfA + l 



(15) 



where Z goes from 1 to iVf. \l) corresponds to eigenvalue 
di (which is the same as n^ when Np = Na) and it is in 
subsystem A. \l) corresponds to eigenvalue l — di and it 
is in subsystem B. Now, by sticking these modes together 
we obtain a mode in whole system (Klich eigenmodes): 



Nf 



10 



'\i)^ + ^^i^d,\if = j2ui'^\k) 



(16) 



k=l 



We note in general the Klich eigenmodes are not the 
same as the occupied eigenmodes of Hamiltonian ma- 
trix h in Eq. (l5|). We illustrate this by calculating 
the Np X Np overlap matrix between these modes, in 
which each row represents an eigenmode of matrix h 
and each column represents one particle Klich eigen- 
mode. Table IT] is the overlap matrix of a sample with 
N — 30, Np = 15, a = 0.1 (see below for definition of 
a and its physical meaning). In this case of strong dis- 
order, each Klich eigenmode has large overlap (s) with 
one (or at most two) eigenmode(s) of matrix h, and one 
can (almost) identify a one-to-one correspondence be- 
tween them. Table |ll] is overlap matrix of a sample with 
N = 30, Np = 15,a = 0.9. In this case of weaker disor- 
der, such a correspondence is not possible. 



P{J) = const. J-^+", 



(17) 



with a slowly decreasing with energy scale. We thus use 
this type of distributions of J's for our Hamiltonian in 
Eq. ([2]) or Eq. Q, with smaller a corresponding to 
broader distribution (on logarithmic scale) and stronger 
disorder. Fig. [2] shows an example of the random singlet 
state generated by the RSRG procedure. 

Such approximate ground state is made of a direct 
product of singlet pairs, and only those singlets that cross 
the boundary contribute to entanglement. In the fermion 
language introduced earlier, there is one fermion living on 
the two sites that form a singlet, with equal probability. 
Thus the entanglement energy of the eigenmodes of h^ 
for subsystems A would only take 3 possible values for 
such a state: +00, —00 and 0; the last corresponds to 
singlets formed by spins on opposite sides of the bound- 
ary, with entangled modes having equal weight on the 
two sites of these spins, and zero weight elsewhere. We 
will show that this is not exactly the case for the exact 
ground state. 



IV. RESULTS 

As mentioned before, two important predictions of 
RSRG method are as follows. First, entangled modes 
are ultra localized at entangled spin sites; second, en- 
tanglement energy of these modes are exactly zero, and 



that corresponds to n 



— 0.5. In this section 



we show that these predictions are not accurate; on the 
other hand the entanglement energy does approach zero 



TABLE I: Overlaps between Klich eigenmodes and filled eigenmodes of original Hamiltonian {Nf x Nf matrix) when A^ = 
30, Nf = 15, a — 0.1. Each row represents an eigenmode of matrix h and each column represents one particle Klich eigenmode. 
One can find a one-to-one correspondence between them in most cases. 
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TABLE If: Overlaps between Klich eigenmodes and filled eigenmodes of original Hamiltonian {Nf x Nf matrix) when 
N — 30, Nf ~ 15, Q = 0.9. Each row represents an eigenmode of matrix h and each column represents one particle Klich 
eigenmode. The overlaps spread among quite a few eigenmodes, and a one-to-one correspondence between them is not possible. 
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A. Entangled modes 



FIG. 1: A schematic plot real space renormalization group 
(RSRG) procedure. Each solid line with coupling J connects 
two spins which are depicted by filled circles. The dashed line 
is an effective coupling generated by RSRG. 



asymptotically in the limit of large separation between 
the entangled modes, reflect the asymptotic exactness of 
the RSRG method. 



Particle-hole symmetry of Hamiltonian Eq. ^ forces 
the entanglement energies to be in ± pairs^^J If we choose 
number of spin sites as: A^ = 2x odd integer, then N^ = 
N/2 is odd, we will have at least one pair of entangled 
modes whose entanglement energy is exactly zero. In this 
subsection we focus on these modes. 

Fig. [3] shows a sample with TV = 102, and we choose 
Np = Na = Nb = N/2 = 51. J's are generated using 
the distribution function of Eq. ( 17 1 with a = 0.1 (strong 



disorder). Panel (a) shows the singlet pair formation ac- 
cording to RSRG procedure. We can see there is just one 
singlet pair that crosses the boundary, and contributes to 
entanglement. On the other hand, we can also calculate 
eigenmodes and eigenvalues of entanglement Hamiltonian 
exactly. Panel (b) shows the zero (entanglement) energy 
mode, the Klich eigenmode with n'^ = n^ = 0.5. We 
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FIG. 2; [Color online] An example of random singlet ground 
state formed in a spin chain oi N = 60, with blue curves con- 
necting spins forming singlet pairs. Only those singlets that 
cross the boundary (red vertical line) contribute to entangle- 
ment. 



see that this Klich eigenmode is indeed strongly local- 
ized around the two spins that are entangled (consistent 
with RSRG). 

For weaker disorder RSRG method may not be as re- 
liable. Fig. [4] shows a similar sample of Fig. [3] but 
with a = 0.9. The RSRG method generates almost the 
same singlet pair configuration [see panel (a)]. However 
in panel (b) we see the zero (entanglement) energy Klich 
eigenmode is much more spread out compared to the pre- 
vious case, even though it still peaks at the location of 
the two spins RSRG predicts to form the entangled pair. 
In such cases (and including the earlier strong disorder 
example) we should view the entanglement as coming 
from two effective spins that actually spread over some 
finite spatial extent, whose profile are described by these 
entangled modes. This is appropriate as the profiles of 
these entangled modes are not sensitive to the position 
of the boundaries of these systems (see panel c of Figs. Is] 
andEl), indicating that they are intrinsic properties of the 
degrees of freedom that are entangled over long distance. 

To illustrate the point above quantitatively, we con- 
sider the same samples of iV = 102, but choose different 
locations of the boundary: Na = 19, 39, 51, 69, 79, 89, 
and calculate the overlaps of zero energy Klich eigen- 
modes among them, forming a 6 x 6 overlap matrix. The 
results are presented in Table [III] for a = 0.1 and in Table 
|IV| for a = 0.9. RSRG predicts the same entangled spin 
pair for them, except for the first choice of Na — 19. 
As a result the overlaps are all close to 1 (closer to 1 
for the stronger disorder case a = 0.1), except for those 
involving Na = 19. 

To quantify the spatial extent of such effective spins, 
we use inverse participation ratio (IPR) of the corre- 
sponding Klich eigenmode: 



FIG. 3: [Color online] (a) Singlet bond configuration gener- 
ated by RSRG method for an A'^ = 102 sample with a = 0.1 
(strong disorder). Panels (b) and (c) are zero energy Klich 
eigenmodes corresponding to two different boundaries. Panel 
(b): Zero (entanglement) energy Klich eigenmode which is 
strongly localized at the entangled spins (consistent with 
RSRG prediction), with Na = 51 (red line), and panel (c) 
corresponds to Na = 69 (green line). 



IPR 



EJV'. 



(18) 



For example, if a wave function is localized over just 
one site, ^pi = dij, then IPR= 1, and if a wave function 
spreads equally over all sits, ■0, = 1/\/N, then IPR= N; 
thus IPR is a measure of how much the wave function 
spreads. In particular, if RSRG were exact, it would 
predict IPR=2 for the Klich modes that correspond to 
the singlet pair contributing to entanglement. It is thus 
reasonable to identify ^^^ as the spatial extent of an 
effective spin that contribute to entanglement, through 
formation of a singlet with another effective spin in the 
other subsystem. 

To see how IPR depends on a, we fix the chain length 
to be A?^ = 2 X 41 = 82, Na = Nb = Np = 41, and vary 
a between 0.1 and 1. We only keep samples in which we 
have only one singlet bond crossing the boundary (and 
contributing to entanglement) according to RSRG. Re- 
sult is shown in Fig. [S] We see IPR is very close to 2 
for strong disorder, indicating the entangled modes are 
strongly localized on individual spins (as RSRG predicts). 
On the other hand it is also clear that IPR increases with 
increasing a (or decreasing disorder strength), indicat- 
ing the entangled modes are spreading out over multiple 
spins, a piece of physics RSRG is unable to capture. 

We find that the spatial spreading of the entangled 
modes is determined by the disorder strength, and in- 
sensitive to system size or length of the entangled bond 
(which typically grows with system size). To demonstrate 



TABLE III: Overlaps between zero energy Klich eigenmodes for a sample with q = 0.1, A^ = 102, A^f = 51 corresponding 
to different choices of boundary, Na ~ 19,39,51,69,79,89. The first one, Na = 19, corresponds to a different entanglement 
configuration. 
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TABLE IV: Overlaps between zero energy Klich eigenmodes for a sample with a = 
to different choices of boundary, Na = 19,39,51,69,79,89. The first one, Na = 19, 
configuration. 
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FIG. 4: [Color online] (a) Singlet bond configuration gener- 
ated by RSRG method for a sample with a = 0.9 (weaker 
disorder). Panels (b) and (c) are zero energy Klich eigen- 
modes corresponding to two different boundaries. Panel (b): 
Zero (entanglement) energy Klich eigenmode which is peaked 
at the entangled spins predicted by RSRG but much more 
spread out, A^a = 51 (red line) and panel (c) corresponds to 
Na = 69 (green line). 



FIG. 5: [Color online] Inverse participation ratio (IPR) of the 
zero (entanglement) energy Klich modes versus a, averaged 
over 1000 samples. The chain length is fixed to be 82. Stan- 
dard error is the red line. 



B. Entanglement energy 



this point we calculate the average IPR of the zero (en- 
tanglement) energy Klich mode versus number of spins, 
N, for fixed a. Again we only keep samples with only one 
singlet bond crossing the boundary according to RSRG. 
result is shown in Fig. |6] For both strong and weak 
disorder, we find little dependence on size N. 



Now, we focus on entanglement energy. If we choose 
the number of spin sites, N, to be 2x even integer, then 
Na = N/2 is even, RSRG would predict we will have an 
even number of entangled pairs of spins (possibly zero) . 
In the following we focus on samples with exactly two of 
such entangled pairs based on RSRG, which then predicts 
there should be two modes whose entanglement energies 
are exactly zero, while all other entanglement energies 
should be ±00. We find that while never literally true. 
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FIG. 6; [Color online] IPR of zero entanglement energy Klich 
modes averaged over 500 samples, for two values of a, versus 
number of sites A^. Green line is average IPR for strong disor- 
der, a — 0.1, and blue line is average IPR for weak disorder, 
a = 0.9. Standard error is the red line. 



FIG. 7: Panel (a): Singlet bond configuration generated by 
RSRG method for an A^ = 100 sample with a = 0.1 (strong 
disorder). Panels (b) and (c): Entangled modes whose en- 
tanglement energies are closest to zero. They are strongly 
localized at the entangled spins (consistent with RSRG pre- 
diction) . 



the above is close to being the case for strong disorder, 
but not quite so for weaker ones. As example we consider 
a system oi N = 100, Np = Na = Nb = N/2 = 50. 
Table [V] lists the entanglement energies of subsystem A, 
e^, and their correspon ding probabilities, nf, when we 

lists e^ and ni} when we use 



VI 



use a = 0.1; while Table 
a = 0.9. We can see from comparison that for strong 
disorder, there is a pair of entanglement energies very 
close to zero, while (most) others have big magnitudes, 
consistent with RSRG prediction. For weaker disorder, 
on the other hand, the lowest entanglement energy is not 
as close to zero, and there are more entanglement energies 
that are of order 1. Fig. [7](for strong disorder) and Fig. 
Is] (for weak disorder) shows RSRG generated singlet pair 
configuration for these two samples, and the two Klich 
eigenmodes with entanglement energy closest to zero. 

To better quantify this point, we calculate smallest 
(in magnitude) entanglement energy averaged over many 
samples, for a range of a's. Fig. [o] shows the (averaged) 
smallest entanglement energy versus a, for samples of 
iV = 60, Np = Na^ Nb = N/2 = 30. As we can see in 
this figure, this energy is close to zero in strong disorder, 
while as we go to weaker disorder it increases, indicating 
the decreasing reliability of RSRG. 

RSRG is expected to be asymptotically exact in the 
limit of long distance and low-energy, at least for ther- 
modynamic properties.'^ Here we would like to test if it 
is also the case for entanglement energy. To this end we 
calculate the smallest entanglement energy versus bond 
length for a specific a. Since we focus here on samples 
with two singlet bonds crossing the boundary, we identify 
the smallest (in magnitude) entanglement energy with 
the bond length of the longer of the two bonds. Since the 
distribution of the entanglement energies is very broad on 
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FIG. 8: Panel (a): Singlet bond configuration generated by 
RSRG method for a.n N = 100 sample with a — 0.9 (weak 
disorder). Panels (b) and (c): Entangled modes whose entan- 
glement energies are closest to zero. They are peaked at the 
entangled spins predicted by RSRG but much more spread 
out. 



logarithmic scale for large size systems, in Fig. 10 we plot 
the average of log e vs. logarithm of bond length. The 
range of the latter is divided into 40 segments, with log e 
within the same segment averaged over. We have also 
included data from samples with different systems sizes 
{N) in the averaging, as we find e depends on bond length 
only. We find, qualitatively, the entanglement energy de- 
creases with increasing bond length, and approaches zero 



TABLE V: Some of e^. 's and corresponding n^ 's for the case of strong disorder (q = 0.1) and TV = 100, A'^f ~ Na ~ Nb 
N/2 = 50. 
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TABLE VL Some of ej, 's and corresponding n^ 's for the case of weak disorder (a = 0.9) and N — 100, Nf ~ Na ~ Nb 

N/2 = 50. 
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FIG. 9: [Color online] Average smallest (in magnitude) en- 
tanglement energy e versus a for a sample oi N = QQ, Nf = 
Na = Nb = N/2 = 30. Average is calculated over 200 sam- 
ples. Red line is standard error. For strong disorder e is close 
to zero and as we approach weak disorder regime it goes up. 



as bond length — > oo. This is consistent with asymptoti- 
cal exactness of RSRG. More quantitatively, we find that 
for two different disorder strengths, a — 0.9 and a — 0.6, 
< log e > depends on logarithm of bond length linearly 
(beyond certain length scale), indicating 



FIG. 10: [Color online] Average logarithm of smallest en- 
tanglement energy (e) versus logarithm of bond length. For 
a — 0.9 we study samples ranging from size A'" = 100 to 
A^ = 1000 with interval of AA = 100. For a = 0.6, we study 
samples with A^ = 100, 200 and A^ = 300. For each TV we cal- 
culate smallest (in magnitude) e and the corresponding bond 
length for 500 samples. Logarithm of bond length is divided 
into 40 segments, with log e within each segment averaged 
over. The blue line is best linear fit for a = 0.9 data with 
slop of —8.3 with standard diviation 0.5. Green line is the 
best linear fit for the case of a = 0.6 data with slop of —8.5 
with standard deviation 0.6. 



SUMMARY AND CONCLUDING REMARKS 



^typical 



L- 



a = 8.3 ±0.5, a = 0.9, 
a = 8.5 ±0.6, a = 0.6, 



(19) 



where L is bond length, etypicai — e^'°S'^^ is the typi- 
cal value of entanglement energy (in contrast to average 
value, which is often dominated by rare fluctuations for 
broad distributions). Such power law behavior is con- 
sistent with the quantum critical nature of the random 
singlet phase. The corresponding exponent is the same 
for the two cases within error bars, indicating its univer- 
sality. 



In this paper we have studied the entanglement spec- 
trum and in particular, the pairs of entangled modes of 
random spin- 1/2 XX chains using a free fermion map- 
ping, and compared them with predictions of real space 
renormalization group (RSRG) treatment. We find that 
RSRG results are qualitatively valid, but not quantita- 
tively accurate, especially for modest disorder strength. 
On the other hand, in the large distance limit its pre- 
diction about entanglement energies becomes asymptoti- 
cally exact. It would be interesting in the future to study 
if such asymptotic exactness still holds in the presence of 
relevant perturbations that drives the RSRG flow away 



9 



from the random singlet fixed point, like dimerizationflSI 
On the other hand, the spreading of the entangled modes 
does not decrease with increasing system size; we inter- 
pret the profiles of these modes as the images of the 
effective spins that form singlets in the random singlet 
phase. This demonstrates the usefulness of studying the 
entangled modes, as such effective spins cannot be stud- 
ies quantitatively within RSRG. 

In a broader context, we argue that eigenstates of the 
reduced density matrices contain much useful informa- 
tion about entanglement, just like their spectra which 



are widely studied now. While the example presented 
here is a very simple one, we hope it serves as a starting 
point for future studies of entanglement taking advantage 
of eigenstates of the reduced density matrices. 
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